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• Abstract. IRAM 30m ^^CO(l-O) and ^^C0(2-l) HERA observations are presented for the ram-pressure stripped 

Virgo spiral galaxy NGC 4522. The CO emission is detected in the galactic disk and the extraplanar gas. The 
Q , extraplanar CO emission follows the morphology of the atomic gas closely but is less extended. The CO maxima 

' do not appear to correspond to regions where there is peak massive star formation as probed by Ha emission. 

' The presence of molecular gas is a necessary but not sufficient condition for star formation. Compared to the 

, disk gas, the molecular fraction of the extraplanar gas is 30 % lower and the star formation efficiency of the 

extraplanar gas is about 3 times lower. The comparison with an existing dynamical model extended by a recipe 
for distinguishing between atomic and molecular gas shows that a significant part of the gas is stripped in the 
^ . form of overdense arm-like structures. It is argued that the molecular fraction depends on the square root of the 

QQ ' total large-scale density. Based on the combination of the CO/Ha and an analytical model, the total gas density 

, is estimated to be about 4 times lower than that of the galactic disk. Molecules and stars form within this dense 

gas according to the same laws as in the galactic disk, i.e. they mainly depend on the total large-scale gas density. 
lO ' Star formation proceeds where the local large-scale gas density is highest. Given the complex 3D morphology this 

, does not correspond to the peaks in the surface density. In the absence of a confining gravitational potential, the 

' stripped gas arms will most probably disperse; i.e. the density of the gas will decrease and star formation will 

oo ■ 

o 
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1. Introduction optical radius, and a significant amount of extraplanar gas 

to the west of the galactic disk. The one-sided extrapla- 
The best-studied case of active ram-pressure stripping of ^^^r atomic gas distribution shows high column densities, 
a cluster galaxy is NGC 4522 located in the Virgo cluster, comparable to those of the adjacent galactic disk. The 
Only m this cluster is the resolution of radio telescopes g ^m polarized radio continuum emission shows a maxi- 
sufficient 20" = 1.6 kpcHI) for a detailed analysis of the ^^um at the eastern edge of the galactic disk, on the oppo- 
gas morphology and kinematics. NGC 4522 is a rather gj^g g^^g extraplanar gas and star formation. Since 
small (D25 = 4' = 20 kpc) edge-on Sc galaxy with a ro- ^^le stellar disk is symmetric and undisturbed (Kenney & 
tation velocity of ^ 100 kms \ It is strongly Hi deficient Koopmann 1999), a tidal interaction is excluded as the 
{DEF - 0.6; Helou et al. 1984). Its projected distance ^^^^^^ ^^^^le pecuUar gas distribution of NGC 4522. Thus, 
to the cluster center (M 87) is large (~ 1 Mpc), and its ^^lis galaxy undergoes ram-pressure stripping due to the 
radial velocity with respect to the Virgo cluster mean is galaxy's rapid motion within the hot and tenuous intra- 
high (1150 kms ). Hi and Ha observations (Kenney et (.j^gter gas (ICM) of the Virgo cluster, 
al. 2004; Kenney & Koopmann 1999) show a heavily trun- 
cated gas disk at a radius of 3 kpc, which is ~ 40% of the VoUmer et al. (2006) made a dynamical model that 

includes the effects of ram pressure for NGC 4522. The 

Send offprint requests to: B. VoUmer, e-mail: model successfully reproduces the large-scale gas distribu- 

bvollmer@astro.u-strasbg.fr tion and the velocity field. By assuming a Gaussian distri- 

* Based on IRAM 30m HERA observations bution of relativistic electrons, they obtained the distribu- 

^ We use a distance of 17 Mpc for the Virgo cluster. tion of polarized radio continuum emission, which repro- 
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duces the VLA observations of polarized radio continuum 
emission at 6 cm. The observed maximum of the polar- 
ized radio continuum emission is successfully reproduced. 
The eastern ridge of polarized radio continuum emission is 
therefore due to ram pressure compression of the interstel- 
lar medium (ISM) and its magnetic field. The dynamical 
model and the analysis of the stellar populations of the 
outer gas-free disk using optical spectra and UV photome- 
try (Crowl & Kenney 2006) indicate that the ram pressure 
maximum occurred only ~ 50-100 Myr ago. This scenario 
has one important caveat: the large projected distance 
of NGC 4522 (1 Mpc) to the center of the Vhgo clus- 
ter (M87). Assuming a static smooth ICM and standard 
values for the ICM density and the galaxy velocity, the 
ram pressure at that location seems to be too low by an 
order of magnitude to produce the observed truncation of 
the gas disk. A natural explanation for the enhanced ram 
pressure efficiency is that the intracluster medium is not 
static but moving due to the infall of the M49 group of 
galaxies from behind (Kenney et al. 2004, VoUmer et al. 
2004, VoUmer et al. 2006). In this case the galaxy has just 
passed the region of highest intracluster medium velocity. 

While we know from the Hq observations (Kenney et 
al. 2004, Kenney & Koopmann 1999) that stars are form- 
ing in the extraplanar gas, we do not know the distri- 
bution of molecular gas in these regions. How does the 
complex multiphase interstellar medium respond to ram 
pressure stripping? Can we model the molecular gas con- 
tent during the interaction using simplified recipes? Can 
dense molecular gas decouple from the ram pressure wind 
as suggested for NGC 4438 (VoUmer et al. 2005)? In this 
article we present IRAM 30m 12CO(1-0) and 1^00(2-1) 
HERA observations of NGC 4522 to investigate the fate 
of the stripped gas. 

We present our CO observations in Sec. [2] followed by 
the observational results in Sec. [H The detection of ram 
pressure wind decoupled molecular clouds is reported in 
Sec. 131 In Sec. [5] we compare our CO observations to exist- 
ing Hi and Ha emission distributions (Kenney et al. 2004) 
and to the dynamical model of VoUmer et al. (2006). The 
molecular fraction and star formation efficiencies are dis- 
cussed in Sec. |6]and we give our conclusions in Sec. [71 

2. Observations 

The observations of the CO(l-O) and C0(2-l) lines, with 
rest frequencies of 115.271204 and 230.53799 GHz respec- 
tively, were carried out at the 30 meter millimeter-wave 
telescope on Pico Veleta (Spain) run by the Institut de 
RadioAstronomie Millimetrique (IRAM). The C0(2-l) 
observations used the HERA multi-beam array, with 3 x 
3 dual-polarization receivers, and the WILMA autocorre- 
lator backend with 2MHz spectral resolution. The C0(1- 
0) observations used the single-pixel "AB" receivers and 
the IMHz filterbanks as backends. The spectral resolu- 
tion is 2.6 kms~^ in both cases. The HERA observations 
were made in February and March 2006 and the CO (1-0) 
in November 2006. In both cases, a nutating secondary 




Fig. 1. IRAM 30m pointings on a DSS B band im- 
age. Triangles: ^^C0(2-l) HERA observations. Small cir- 
cles: ^^CO(l-O) observations. Large circles: Binned HERA 
pointings (see Fig. [7]). The large boxes correspond to Fig.[3l 
to Fig. El 



("wobbler") was used with a throw of 180-200 arcseconds 
in order to be clear of any emission from the galaxy. The 
positions observed in each line are indicated in Figs. [Hand 
[lias triangles for C0(2-l) and chcles for CO(l-O). 

Data reduction was straightforward, eliminating any 
obviously bad channels and excluding the spectra taken 
under particularly poor conditions (system temperature 
over 1000 K). Spectra were then summed position by posi- 
tion. System temperatures of the final spectra ranged from 
200 to 500 K on the Ta* scale. All spectra are presented 
on the main beam temperature scale, assuming telescope 
main-beam and forward efficiencies of ?7mb — 0.54 and 
T^fo^ = 0.90 for HERA and Tymb = 0.74 and T^for = 0.95 
for the CO(l-O) line. The spectra near map edges with 
noise levels greater than 28 mK {Tmb scale) are left out 
of Fig.[l At the assumed distance of NGC 4522, 17 Mpc, 
the C0(2-l) and CO (1-0) beams correspond to 0.9 and 
1.7 kpc respectively. In order to convert CO integrated 
intensities into molecular gas masses, we have assumed 
a iV(H2)//co(2-i) ratio of 2 x 10^° H2 mol cm-^ per K 
kms^^. Our conclusions, however, do not depend strongly 
on the -/V(H2)//co ratio within reasonable variations. 
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Fig. 2. 12C0(2-1) HERA pointings (see Fig. [T]) with low 
noise levels on the Hi emission distribution (from Kenney 
et al. 2004). The large boxes correspond to Fig.[3]to Fig. [HI 

3. Disk and extraplanar molecular gas 

3.1. CO spectra 

Fig.lMlshow the ^200(2-1) HERA spectra (resolution of 
11") together with the Hi spectra of Kenney et al. (2004) 
(resolution of 20") and the model Hi and CO emission con- 
volved to the observational resolutions (see Sect. [5]). We 
observe that the Hi emission is more extended than the 
CO emission. Where detected, the velocity of the CO lines 
are close to the Hi velocities. In general, CO linewidths 
are comparable to the Hi linewidths in the galactic disk 
(Fig. HI, but smaller in the extraplanar regions (Fig. \3\ 
[6]). A double-hne profile is observed in the CO and Hi 
line west of the galaxy center at offsets (—23, 0), (—23, 8) 
(Fig. [5]). This kind of line- profile has also been observed 
in NGC 4438 (VoUmer et al. 2005) another Virgo spiral 
galaxy which undergoes ram pressure stripping together 
with a tidal interaction. In the southwestern part of the 
extraplanar gas the Hi profiles show a blueshifted wing 
(Fig. [S]), corresponding to the most strongly pushed gas. 
Whereas the CO peak in this region is aligned with the Hi 
peak, the blueshifted wing is absent in CO. This might be 
partly due to the smaller S/N ratio of the CO data com- 
pared to the Hi data. Even if there is H2 associated with 
the blueshifted wing, we can conclude that the molecular 
fraction in this extraplanar gas is lower in the blueshifted 
wing than in the main line. 

Fig. [7] shows i2C0(2-l) HERA spectra of selected re- 
gions (large circles in Fig. [1] [2]). To obtain a better S/N 




Fig. 3. IRAM 30m HERA 1^00(2-1) spectra of the north- 
eastern box of Fig. [T] (black lines). Hi spectra from Kenney 
et al. (2004) (red lines). CO (green lines) and Hi (blue 
lines) spectra from the dynamical model of VoUmer et al. 
(2006). The velocity scale is from 2150 to 2550 kms'^ and 
the main-beam temperature scale from -40 mK to 80 mK 
(T„jfc scale). 

ratio, the HERA spectra within each region were averaged. 
These selected regions show the following characteristics: 

— CO emission is present in the northwestern extraplanar 
gas at offset (15", 64"). The line is redshifted by - 
20 kms~^ with respect to the Hi line (Fig.[7]left panel). 

— The double profile observed in the Hi line at offset 
(-27", 0) is also present in the CO hue. The CO and 
Hi peak amplitudes and their separation are the same 
(Fig. [7] middle panel) . 

— In the southwestern extraplanar gas at offset 
(-46", -21") the blueshifted wing of the Hi hue pro- 
file has no counterpart in the CO data (Fig. [7] right 
panel). The CO linewidth is significantly smaller than 
the linewidth of the main Hi line. This confirms the low 
molecular fraction of the blueshifted diffuse atomic gas. 
This blueshifted component corresponds to low surface 
density atomic gas discussed in Kenney et al. (2004; 
see their Fig. 10). The lower molecular fraction of this 
gas confirms the claim of VoUmer et al. (2006) that 
this gas has low densities making it more vulnerable 
to ram pressure stripping. 

The ^^CO(l-O) spectra (resolution: 21") are shown to- 
gether with the convolved CO (2-1) HERA spectra in 
Fig.il For aU but offset (20", 0), the CO(l-O) line closely 
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Fig. 4. IRAM 30m HERA CO (2-1) spectra of the south- 
eastern box of Fig. [T] (black hnes). Hi spectra from Kenney 
et al. (2004) (red lines). CO (green lines) and Hi (blue 
lines) spectra from the dynamical model of Vollmer et al. 
(2006). The velocity scale is from 2150 to 2550 kms'^ and 
the main-beam temperature scale from -40 mK to 80 mK 
{Tmb scale). 




Fig. 5. IRAM 30m HERA 1^00(2-1) spectra of the north- 
western box of Fig. [T] (black lines). Hi spectra from Kenney 
et al. (2004) (red lines). CO (green lines) and Hi (blue 
lines) spectra from the dynamical model of Vollmer et al. 
(2006). The velocity scale is from 2150 to 2550 kms^^ and 
the main-beam temperature scale from -40 mK to 80 mK 
{Trab scale). 



follows the C0(2-l) line with a line ratio of 0.7-0.9. Thus, 
the density and temperature of the molecular gas in the 
western extraplanar regions is probably not significantly 
different from those of the molecular gas in the galac- 
tic disk. At the eastern edge we observe a lower C0(2- 
1)/CO(1-0) ratio (~ 0.5). Given the sharpness of the gas 
distribution at the eastern edge and the pointing uncer- 
tainty of the telescope at 115 GHz, this small line ratio 
could be due to an offset of the CO (1-0) pointing to the 
west, i.e. closer to the major axis, with respect to the 
C0(2-l) pointing. 

4. Wind-decoupled molecular gas 

In another Virgo spiral galaxy, NGC 4438, Vollmer et 
al. (2005) found CO emission not associated with any 
Hi emission. NGC 4438 underwent a tidal interaction 
~ 100 Myr ago (Combes et al. 1988) and now under- 
goes severe ram pressure stripping (Vollmer et al. 2005). 
A narrow CO line was detected in the northern tidal arm 
of NGC 4438, with apparently no associated Hi. Since the 



velocity of the CO line corresponds to that of the stellar 
component (determined using a dynamical model of the 
tidal interaction), Vollmer et al. (2005) claimed that these 
molecular clouds were too dense to be affected by ram 
pressure, i.e. that they decoupled from the ram pressure 
wind. 

Based on these findings we searched for CO lines in 
regions devoid of any Hi emission (i) in the ram pres- 
sure stripped outer galactic disk (Fig. [9]) and (ii) between 
the northern galactic disk and the extraplanar Hi emission 
(Fig. [T0|) . We only detect one CO line in the northern part 
of the galactic disk. The most prominent CO line is de- 
tected at the position (30, 68); see the little box in Figs.[l] 
and[2|), i.e. between the stellar disk and the stripped ex- 
traplanar atomic gas. The radial velocity of the CO lines 
is ^ 2250 kms^^, close to the velocity of the stellar com- 
ponent. The Hi line of the stripped extraplanar Hi (Fig. [3]) 
is blueshifted by ~ 30 kms~^. We observe a clear absence 
of a CO line in the gas-free southwestern part of the disk 
(Fig. [9] right panel), for which an explanation is proposed 
in Sect. 6.1. 
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Fig. 6. IRAM 30m HERA 12C0(2-1) spectra of the south- 
western box of Fig. [T] (black hnes). Hi spectra from Kenney 
et al. (2004) (red lines). CO (green hnes) and Hi (blue 
lines) spectra from the dynamical model of VoUmer et al. 
(2006). The velocity scale is from 2150 to 2550 kms'^ and 
the main-beam temperature scale from -40 mK to 80 mK 
{Trnb scale). 




Fig. 7. Binned IRAM 30m HERA 12C0(2-1) spectra (see 
large circles on Fig. [T]). Hi spectra from Kenney et al. 
(2004) (red hnes). CO (green hnes) and Hi (blue hnes) 
spectra from the dynamical model of VoUmer et al. (2006). 
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Fig. 9. Average ^^C0(2-l) spectra of the northeastern 
(45, 68) and the southwestern (—38, —61) ends of the 
galactic disk, where the atomic hydrogen has been re- 
moved by ram pressure. No Hi emission is detected nor 
does the model predict any gas. The spectra are averaged 
over a region of 20". 
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Fig. 10. Single ^^00(2-1) HERA spectrum in the north- 
east that is located outside the Hi distribution. The posi- 
tion is marked with an additional box in Fig. [2l 



5. Comparison with the dynamical model 

In this section we will compare the molecular gas and 
Ha distribution to the dynamical model of VoUmer et al. 
(2006). 
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Fig. 8. Black: ^^CO(l-O) spectra (see small circles on 
Fig. [J). Red: convolved i2C0(2-l) HERA spectra. 



5.1. Molecular gas 

In order to estimate the H2 column density distribu- 
tion, via a zero moment (integrated intensity) map of the 
C0(2-l) emission, we proceed as below. For the positions 
where the S/N of the Hi is higher than 2.5(7, the linewidths 
of the Hi spectra of these positions are determined and the 
CO spectra integrated over the Hi velocity range. This 
leads to the CO emission distribution map of Fig. [11] (top 
panel) . 

To compare the Hi and the CO emission distributions 
with the dynamical model of VoUmer et al. (2006), we 
assume that the molecular fraction of the gas depends on 
the local gas density. In a first simple approach we assume 
that the gas is bimodal, i.e. entirely molecular at densities 
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Fig. 11. Top panel: ^^C0(2-l) emission distribution (solid 
line) on the Hi emission distribution (dotted line). The 
plus signs show the IRAM HERA pointings which were 
used to construct the CO moment map. Middle 
panel: model CO distribution (solid line) on the model 
Hi emission distribution (dotted line). The gas is as- 
sumed to be entirely molecular for gas densities p > 
0.03 Mqpc"'^. Lower panel: model CO distribution (solid 
line) on the model Hi emission distribution (dotted line). 
The molecular gas fraction is calculated using /moi = 
a//9/(0.5 Mopc-3). The letters refer to the regions shown 
in the deprojection in Fig. [T31 



> 0.03 M0pc~'^ = 1 cm~3 and entirely atomic at lower 
gas densities. The resulting model gas distribution maps 
were convolved to the observational resolutions (middle 
panel of Fig. [TT|) . In a second approach we assume that the 
molecular fraction depends linearly on the square root of 
the gas density /,„oi = v^p/(0.5 Mopc"^). Moreover, the 
molecular fraction cannot exceed unity. In Sec. 16. 21 we will 
give a motivation for this dependency. The resulting model 
CO emission distribution is shown in the lower panel of 
Fig.im 

The comparison between the observed and the simu- 
lated CO emission distribution shows the following simi- 
larities: 

— The galactic disk is the most prominent feature. There 
is more molecular gas of higher surface densities in the 
outer part of the gas disk in the southwest than in the 
northeast. 

— There is CO emission in the extraplanar regions. 

— The model using a molecular fraction proportional to 
the square root of the gas density better reproduces the 
CO morphology of the extraplanar regions: the north- 
ern and southwestern CO emission regions are centered 
on the peaks of the Hi emission and there is a spatially 
separated CO arm between the southwestern and the 
disk CO emission. 

On the other hand, we observe the following disagreement 
between the model and our observations: 

— The model CO emission distribution shows a central 
hole. This is due to the initial conditions which had an 
initial gas hole for computational reasons. 

— The northern part of the observed CO disk shows 
emission of lower surface brightness. The Hi emission 
(Kenney et al. 2004) shows the opposite trend: high 
column density gas is found to the north. 

— The model CO emission to the west of the galaxy cen- 
ter is more extended than it is observed. This is also 
the case for the Hi emission (see Sect. [6]). 

We thus conclude that the model reproduces qualita- 
tively our CO observations. The model using a density- 
dependent molecular fraction reproduces the observations 
better than a bimodal molecular fraction. 

The direct comparison between the Hi and CO model 
and observed spectra (Fig. [3] to Fig. [7]) shows good agree- 
ment for high intensities, i.e. high gas densities. However, 
the observed CO and Hi double lines and blueshifted wings 
of the Hi lines in the western extraplanar regions are not 
reproduced by the model. This is due to the constant col- 
umn density of the model gas in the regions affected by 
ram pressure (see Sect.[6|). 

5.2. Star formation 

In typical spiral galaxies the star formation rate fol- 
lows more closely the molecular gas distribution than the 
atomic gas distribution (Wong & Blitz 2002). The Hi sur- 
face density saturates at a value of ^ 10 Mqpc~^ or even 
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Fig. 12. Ha emission distribution (Kenney et al. 
2004; contours) on the ^^C0(2-l) emission distribution 
(greyscale) . 

declines for high star formation rates per unit surface. 
In an unperturbed galactic disk the ISM is confined in 
the gravitational potential of the disk. The ISM is turbu- 
lent and this turbulence is most probably maintained by 
the energy input from SN explosions (see, e.g. MacLow 
& Klessen 2004 or Vollmer & Beckert 2003). Without a 
constant energy supply, turbulence is damped within a 
few Myr (Stone et al. 1998, MacLow 1999). Since the ex- 
traplanar gas of NGC 4522 is no longer confined to the 
potential of the galactic disk, it represents an ideal lab- 
oratory to test if the gravitational potential plays a role 
for the correlation between star formation rate and the 
available molecular/ atomic gas mass. 

To do so, we first present the Ha emission distribu- 
tion overlaid onto the CO emission distribution (Fig. [T2|) . 
There is no major site of massive star formation with- 
out associated CO emission. At the edges and outside the 
galactic disk star formation does not coincide with the 
maxima of CO emission. A high column density of molec- 
ular gas is thus not sufficient to form massive stars. In 
other words, on kpc scales the star formation rate does 
not directly depend on the molecular gas surface den- 
sity. A number of Hii regions are detected at the outer 
edges of the CO emission distribution, at offsets (20", 45"), 
(-10", 20"), (-45", 0), (-45", -30"). 

In the model we assume that the star formation rate 
is proportional to the number of collisions between the 
gas clouds. Numerically, the star formation rate thus de- 
pends on the local number density of the clouds, their cross 
section, and their local 3D velocity dispersion. Since the 
model clouds have a constant surface density (Vollmer et 
al. 2006), their cross sections vary with the cloud mass in 
the following way: Trr^j = Md/S, where rd is the radius. 



Mc\ the cloud mass and S the gas surface density. The 
cloud mass distribution is a power law with an index of 
— 1.5. For an isolated unperturbed spiral galaxy this pre- 
scription leads to a Schmidt law of the form oc S^ '', 
where is the star formation rate per area. For the con- 
struction of a star formation distribution map we store all 
cloud-cloud collisions during 20 Myr before the present 
state of the galaxy. This is twice the timescale for Ha 
emission, chosen to give more collisions and thus better 
statistics. The distribution was then convolved to 0.6 times 
the resolution of the C0(2-l) map. This model star for- 
mation distribution is presented in Fig. [T3] together with 
the Ha emission distribution. We note the following points 
for comparison between the model and observations: 

— As in the observations, our model has two prominent 
Hii regions at the outer edges of the galactic disk, 
the most prominent being the southwestern Hii region. 
These regions are separated by a local minimum from 
the rest of the galactic disk. 

— We observe extraplanar model star formation regions 
close to the disk at the edges and in the middle of the 
galactic disk (see Sect. E]). 

— There is isolated star formation over the whole extra- 
planar Hi emission. In the model the small patches are 
due to multiple collisions of a single massive cloud. 

We thus conclude that the model star formation distri- 
bution qualitatively reproduces the observed Ha emission 
distribution. 

6. Discussion 

6.1. Deprojecting the model 

The model gives us the unique opportunity to deproject 
the gas distribution. One has to keep in mind that the 
extraplanar gas is no longer located within the galactic 
disk, but is a fully 3D feature. In Fig. [TH we present a 
deprojected face-on view of the model gas distribution of 
Fig. [Tl] (lower panel). The corresponding regions are la- 
beled with capital letters. Region A is the tip of a spiral 
arm close to the main gas disk. Region B is a large over- 
density within the stripped gas. Region C represents the 
tip of the most prominent gas (spiral) arm which is mainly 
made of stripped material. Region D is a secondary gas 
arm and region E is an overdensity in the windward, low- 
density side of the prominent gas (spiral) arm. As stated 
above, the gas outside the main gas disk (i? < 30") has 
a fully 3D structure. As can be seen in the lower panel of 
Fig.HH the upper edge of the deprojected gas distribution 
becomes the northeastern edge of the model distribution 
in the observed projection. The prominent gas (spiral) arm 
(D, C) runs in the observed projection vertically from the 
southwestern edge of the galactic gas disk to the end of 
the gas distribution at (40,40) (Fig. [H]). The southwest- 
ern extraplanar Hi emission region is made of relatively 
diffuse gas from the low-density side of the prominent gas 
(spiral) arm. This low-density side has one overdensity (E) 



8 



VoUmer et al.: Stripped molecular gas in NGC 4522 



] r 
NGC 4522 




RIGHT ASCENSION (J2D00) °f'se' (orcsec) 



Fig. 13. Left panel: Hi emission distribution on Ha emission distribution (Kenney et al. 2004). Right panel: model 
Hi distribution (VoUmer et al. 2006; contours) on the model massive star formation distribution (greyscale). Darker 
regions correspond to less massive star formation. 
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Fig. 14. Deprojected model CO distribution (sohd line) 
on the deprojected Hi distribution (dotted line). The let- 
ters correspond to the characteristic regions of Fig. [11] 
The galaxy is moving towards the lower left corner and 
rotates clockwise. 



where star formation can proceed. This is consistent with 
the observed Hi and Ha emission distribution, but we do 
not claim that this is the necessary configuration. Most 
interestingly, we can identify the upbending arm (south 
of D, C on the lower panel of Fig. [TT|) with a chain of 
observed Hii regions (between (12/i33m38s, OOdlOmOOs) 



and (12/i33to37s, 09dlOTO30s) on the left panel of Fig. [13] 
and the corresponding molecular arm on Fig. 1121 We think 
that this is a robust feature. As already stated in Sect. 15. ll 
the outer part of the upbending arm (DEC offset 0" to 
40" in the lower panel of Fig. [TT]) is not present in the 
observations. In deprojection this part of the arm is lo- 
cated between DEC offsets of 20" and 80" (Fig.[14l). Due 
to its low total gas surface density E with respect to the 
inner part (DEC offset —20" to 0") it is more vulnerable 
to ram pressure stripping or evaporation (the acceleration 
a due to ram pressure p is a — p/T,). This is consistent 
with the blueshift and large linewidth of the extraplanar 
low column density Hi which is stripped more efficiently 
than the high column density gas (Fig. 9 of VoUmer et 
al. 2006). Moreover, the CO and Hi spectra of the south- 
ern extraplanar Hi region (Fig. [5]) show that CO is only 
associated with high column density Hi at the highest ve- 
locities. Our numerical model cannot reproduce the more 
efficient stripping of low column density gas, because it 
uses a constant column density for the atomic gas phase 
(see VoUmer et al. 2001 and Fig. 9 of VoUmer et al. 2006). 
We therefore suggest that the outer part of the stripped 
gas arm has been stripped more efRcienctly and now has 
a column density too low to be detected in the Hi obser- 
vations. On the other hand, very dense gas can decouple 
from the ram pressure wind (see Sect. as it is found 
in CO observations of NGC 4438 (Fig. 4 of VoUmer et al. 
2005). 

We conclude that the formation of molecular clouds 
and star formation mainly depend on the large-scale over- 
density of the gas more than on dynamical criteria or the 
overall pressure (see Sect. 6.3). 
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Wind-decoupled molecular gas is only found in the 
northern part of the galactic disk (Sec.|4]). With the help of 
the deprojection (Fig. [T4)) . this can be understood. High 
density gas is stripped from the left border of the gas 
disk. Since the ram pressure wind is rather face-on, the 
gas clouds are pushed to larger heights above the galactic 
plane. At the same time, rotation makes clouds move to- 
wards positions A and B, i.e. the northern part of the disk 
where the wind-decoupled molecular clouds are found. 
This finding suggests that dense giant molecular clouds 
can decouple from the ram pressure wind at early stages 
of the stripping of dense gas from the galactic disk, as 
seems to be the case for NGC 4438 as well. 

6.2. Molecular fraction 

As seen in Sec. lS.ll and Fig. [Til a molecular fraction which 
is proportional to the square root of the total gas density 
reproduces the observed CO emission distribution bet- 
ter than a simple bimodal molecular gas fraction assum- 
ing a gas density cutoff. VoUmer & Beckert (2003) ap- 
proximated the molecular fraction by the ratio between 
the turbulent crossing time scale ^turb = »'ci/wturb and 
the timescale for molecule formation t^oi = ci/Pci, where 
t'turb is the turbulent velocity dispersion within the cloud, 
Pel = p/^v the cloud density, p the total large-scale gas 
density, <i>v the volume filling factor, and a the constant 
of molecule formation a ~ 3 x 10^ yrMopc"'^: 

^turb 



fmol — 



(1) 



^mol Wturb" $V 

The volume filling factor $v is defined by the condition 
that the gas clouds are self-gravitating, i.e. the turbulent 
cloud crossing time ^turb equals the free fall time te of the 
clouds: 



Wturb 



/3 7r$v 
32G/9 



(2) 



where G is the gravitation constant. Inserting Eq. [2]into 
Eq. [T] leads to 



mol 



3tt 
32g' 



(3) 



Thus, the molecular fraction depends on the square root 
of the cloud density pc\ — ^y^p. For the dependence used 
in Sec. IS.ll we assume a constant volume filling factor (see 
Vollmer & Beckert 2003). A density of 0.5 Mopc"^ im- 
plies a volume filling factor of <i>v = 0.03. This is higher 
than the volume filling factors given in Vollmer & Beckert 
(2003). Eq. [T]is a crude approximation which may overes- 
timate the molecular fraction by a factor of 5-10, probably 
because fmoi is substantially overestimated by tturb/tmoi- 
However, we think that its dependence on the physical pa- 
rameters of the gas are valid. Eq. [3] will be used together 
with a similar expression for the star formation efficiency 
in the next section. The observed molecular gas fraction 
decreases from 50 % within the galactic disk to 35 % in the 
extraplanar region (Table [1]) . 



Table 1. Derived masses, molecular fractions, and star 
formation rates/timescales. 





galactic disk 


extrapl. 


total 


frac. extrapl. 


Mm (10" Mq) 


2.5 


1.5 


4 


0.4 


Mh, (10» Mo) 


2.2 


0.8 


3 


0.25 




0.88 


0.53 


0.75 




SFR (Meyr-^) 


0.1 


0.015 




~ 0.14 




2.2 


5.3 






t?or(Gyr) 


4.7 


15.3 







6.3. The efficiency of extraplanar star formation 

Does the star formation efficiency (SEE) change once the 
gas has left the confining gravitational potential of the 
galactic disk? The role of large-scale processes in pro- 
voking star formation is subject to debate, with many 
"recipes" providing reasonably similar fits to observations 
- e.g. the "Toomre" criterion (Kennicutt 1989), a pressure- 
based criterion (Blitz & Rosolowsky 2006), or a basic 
Schmidt (1959) law. The problem is not so much predict- 
ing the behavior in spiral disks but understanding what 
governs large-scale star formation in general in order to be 
able to understand other environments, typically those at 
intermediate and high redshifts. In classical dwarf galax- 
ies, it is difficult to study the SEE because the low metal- 
licity makes the measure of the H2 mass uncertain. In 
Tidal Dwarf Galaxies, morphologically similar but with 
higher metallicities, Braine et al (2001) found that, curi- 
ously, the SFE was not identifiably different from spiral 
galaxies typically 100 times more massive. In the post- 
collision Taffy galaxies or UGC 813/816 system, on the 
other hand, the SFE in the bridge gas is much lower than 
within spiral disks (Braine et al. 2003, 2004). 

The SFE can be defined with respect to the molec- 
ular gas mass available or with respect to the total gas 
mass available, either SFE~^ = t|™ = M(H2)/Msfr or 
SFE"^ = if™ = M(H2 + HI)/MsFR- The observed gas 
masses, molecular fractions and star formation rates and 
timescales are presented in Table [TJ Whereas close to half 
of the Hi is found beyond the galactic disk, this ratio de- 
creases to 1/4 for the molecular gas and to 1/7 for the 
Ha emission. The averaged star formation timescale based 
on the molecular or total gas mass increases from 2.2 or 
4.7 Gyr within the galactic disk to 5.3 or 15.3 Gyr respec- 
tively within the extraplanar region. The star formation 
efficiency thus decreases by a factor of ~ 3 between the 
disk and the extraplanar region. 

In the framework of the model of Vollmer & Beckert 
(2003) the local star formation rate is given by 



P^ 



(4) 



where <I>v is the probability of finding a self-gravitating 
cloud, i.e. the volume filling factor of self-gravitating 
clouds. Inserting the expression for the free fall time of 
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Eq. [5] into Eq. [3] yields the following expression for the 
star formation timescale which corresponds to the inverse 
of the star formation efficiency: 



p* 



37r 
32G 



(5) 



Thus, the star formation timescale depends on the inverse 
of the square root of the total large-scale density and the 
volume filling factor. 

From Table [T] we obtain: 



disk 
mol 



fcxt 
•> mol 



1.3 , and 



t1 



3.3 



leading to 



rT^disk „ 

*v ^2.5, and ^=4.3 



^oxt 



Pcxt 



(6) 



(7) 



We therefore suggest that the observed decrease of the 
star formation efficiency by a factor of 3 in the extra- 
planar region, together with a lower molecular fraction, 
is due to a higher volume filling factor of self-gravitating 
clouds in the galactic disk of NGC 4522. These clouds 
are about twice as dense as their counterparts in the ex- 
traplanar regions. The overall density in the extraplanar 
region is about 4 times lower than that of the galactic 
disk, presumably because the extraplanar gas is no longer 
confined by the gravitational potential of the disk. This 
is supported by the lower CO(2-1)/(1-0) fine ratio in the 
extraplanar gas (Fig. 8). However, the gas is still confined 
by the hot intracluster medium and partially compressed 
by ram pressure. The mixture of ram-pressure and rota- 
tion create zones where the atomic gas is dense enough 
to be gravitationally bound, become molecular, and form 
stars. 

We conclude that the stripped ISM still forms 
molecules and stars in a way not distinguishable from disk 
star formation (Eqs. [T] and [4]) as long as the overall gas 
density is high enough to form bound clouds. The ulti- 
mate fate of the stripped gas is probably ionization and 
evaporation, without star formation for the low-density 
gas, and after a generation of stars, which then disperse 
the remaining dense gas, for the initially denser gas. 

As the simulations show (Fig. [Tl)l part of the gas is 
stripped in relatively dense arms whose mean density is 
about 4 times lower than that of the galactic gas disk. 
Since these gas arm are only confined by the hot intraclus- 
ter medium, they might ultimately disperse giving rise to 
a large low surface density tail as observed in NGC 4388 
(Oosterloo & van Gorkom 2005). 

6.4. Comparison with radio continuum observations 

Recently, Murphy et al. (2008a, 2008b) compared Spitzer 
24 /xm emission with 20 cm radio continuum maps. The 
radio-FIR correlation is used to predict the radio emission 
from the Spitzer 24 /j,m emission. They found a radio defi- 
cient region at the eastern outer edge of NGC 4522 's disk 



where ram pressure is pushing the interstellar medium. 
Since the 24 ^m dust emission is associated with molecu- 
lar gas, we compare in Fig. [15] the 6 cm radio continuum 
emission from Vollmer et al. (2004) with the CO emis- 
sion distribution. The extraplanar 6 cm radio continuum 
emission extends farther to the west than the CO emis- 
sion. The radio emission of the northeastern disk extends 
in the region where we found CO emission without asso- 
ciated Hi emission (Fig.[5|). The northern extraplanar CO 
and Hi emission does not show associated radio continuum 
emission. We convolved the 6 cm radio continuum and CO 
maps to a resolution 20" x 20" and computed a ratio map 
which is shown in the right panel of Fig. [9l As expected, 
the 6 cm/CO ratio is highest in the extraplanar region. 
The smallest 6 cm/CO ratios are found at the extremities 
of the disk. We observe a maximum of the 6 cm/CO at 
the eastern edge of the galactic disk where Murphy et al. 
(2008a/b) detect a radio deficient region. Since we do not 
expect the 24 /xm emission to be proportional to the CO 
emission, the interpretation of the result is difficult. In ad- 
dition, we think that our data are not sensitive enough to 
draw a firm conclusion on the 6 cm/CO ratio map. 



7. Conclusions 

We present IRAM 30m i2cO(2-l) HERA and ^2C0(1- 
0) observations of the ram pressure stripped Virgo spiral 
galaxy NGC 4522. We directly compare the CO spectra 
to the Hi data cube of Kenney et al. (2004). In a sec- 
ond step a CO emission distribution map is produced in 
the regions where Hi is detected. The CO emission distri- 
bution is compared to Ha observations of Kenney et al. 
(2004) and to the model distribution of molecular gas de- 
rived from dynamical simulations of Vollmer et al. (2006). 
A map of the distribution of star formation based on the 
numerical cloud-cloud collisions is produced which is then 
compared to the Ha emission distribution. The 3D model 
snapshot allows us to deproject the observed features and 
understand their origin. From this work we conclude that 

1. CO emission is associated with the extraplanar atomic 
gas. The morphology of the molecular gas closely fol- 
lows but is less extended than the Hi morphology. 

2. In the northern part of the galactic disk we find CO 
emission without an Hi counterpart. We interpret this 
detection as wind-decoupled molecular clouds as ob- 
served in NGC 4438 (Vollmer et al. 2005). 

3. In the extraplanar region CO emission is always asso- 
ciated with sites of massive star formation as probed 
by Ha emission. At the resolution of our observations, 
there is no correlation between the CO and Ha peaks. 

4. A model using a molecular fraction proportional to the 
square root of the gas density qualitatively reproduces 
our CO observations. 

5. The model star formation distribution, which is nu- 
merically based on cloud-cloud collisions, qualitatively 
reproduces the observed Ha emission distribution. 
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Fig. 15. Left panel: Contours of the 6 cm radio continuum emission from VoUmer et al. (2004) on the CO emission 

distribution (greyscale). The contour levels are (2,4,6,8,10,20,30,40,50,60,80,100) x 20 ^Jy/beam. The resolution 
is 15" X 15". Right panel: Ratio between the 6 cm and CO emission. Brighter regions have higher CO/6 cm ratios. 



6. The deprojection of the model extraplanar gas shows 
that a significant part of the gas is stripped in the form 
of relatively dense arms. 

7. The formation of molecular clouds, and subsequent 
star formation, occurs at peaks in the large-scale vol- 
ume density of the gas, with no clear difference with 
respect to the disk despite the very different condi- 
tions (i.e. stellar density dominates in the disk but is 
negligible in the extraplanar material). 

8. In the disk gas, the molecular and atomic fractions 
are about equal whereas in the extraplanar gas, there 
is twice as much HI as H2, assuming a standard 
N{'R2)/Ico conversion ratio. 

9. The star formation efficiency of the extraplanar gas is 
about 3 times lower than that of the galactic disk. 

10. Using the analytical framework of Vollmer & Beckert 
(2003) we find that the overall total gas density and 
volume filling factor of self-gravitating clouds are fac- 
tors of 4 and 2.5 lower respectively compared to the 
galactic disk. 

In the early phases of ram pressure stripping (~ 50 Myr 
after peak ram pressure; Vollmer et al. 2006) a significant 
part of the stripped gas is in the form of relatively dense 
arms. At the same time some very dense molecular clouds, 
representing a tiny fraction of the stripped gas, can decou- 
ple from the ram pressure wind. Molecules and stars form 
within the stripped dense gas according to the same laws 
as in the galactic disk, i.e. they mainly depend on the over- 
all total gas density. Star formation proceeds where the 
local large-scale gas density is highest. Given the complex 
3D morphology this does not necessarily correspond to the 
peaks of the surface density. In the absence of a confining 



gravitational potential these stripped gas arms will most 
probably disperse, i.e. their density will decrease and star 
formation will cease. 
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